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We show that every correctable subsystem for an arbitrary noise operation can be recovered by a 
unitary operation, where the notion of recovery is more relaxed than the notion of correction insofar 
as it does not protect the subsystem from subsequent iterations of the noise. We also demonstrate 
that in the case of unital noise operations one can identify a subset of all correctable subsystems 
— those that can be corrected by a single unitary operation — as the noiseless subsystems for the 
composition of the noise operation with its dual. Using the recently developed structure theory for 
noiseless subsystems, the identification of such unitarily correctable subsystems is reduced to an 
algebraic exercise. 

PACS numbers: 03.67.Pp, 03.67.Hk, 03.67.Lx 



A basic principle in quantum computation, communi- 
cation and cryptography is the requirement that if quan- 
tum information encoded in a physical system is pre- 
served during some storage or transmission process, then 
the information is represented at every time on a sub- 
system of the overall system. This heuristic concept 
has been formalized as part of the "subsystem principle" 
for realizing quantum information EI Q], and provides 
a central motivation for quantum error correction tech- 
niques. Recently a framework for studying correctable 
subsystems for arbitrary quantum operations was intro- 
duced in 0, Q under the moniker "operator quantum 
error correction". The approach includes standard quan- 
tum error correcting (subspace) codes 0, 0, 0, 0, 
and decoherence- free/noiseless subspaces and subsystems 
[3 El IH 13 13 El 111 Although it does not yield 
new codes, this app roach enables more efficient correc- 
tion procedures yjj Eg, E3 ■ 

Intuitively, the subsystem principle implies that if a 
subsystem is correctable for a given noise operation, then 
at the output stage of the transmission there should be 
an identifiable copy of the subsystem within the output 
Hilbert space. Our first main result in this paper formal- 
izes this notion and thereby contributes to the rigorous 
formulation of the subsystem principle for quantum in- 
formation. Specifically, we show that every correctable 
subsystem for an arbitrary noise operation only suffers a 
change in representation due to that operation. A uni- 
tary operation is therefore sufficient to recover from the 
noise, in the sense of restoring the quantum information 
to the subsystem in which it was encoded, but due to the 
fact that complementary subsystems are not necessarily 
preserved by the noise, the system is not necessarily pro- 
tected from subsequent iteration of the noise operation. 
We therefore say that the subsystem is only recovered, 
but not corrected, by a unitary operation. 

Quantum error-correcting codes are typically designed 



with particular noise models in mind. Given the myriad 
different physical implementations for quantum informa- 
tion processing and communication, and the diverse na- 
ture of the noise that may afflict these, techniques for 
identifying codes for arbitrary noise models are clearly 
of interest. To date, this has only been achieved for 
passive quantum error correction. It is accomplished by 
identifying the decoherence-free/noiseless subsystems for 
the noise model. A recently developed structure theory 
for passive error correction 0, [2(J reduces the identifi- 
cation of such codes to an exercise in matrix algebra. 
(The algorithm is polynomial in the dimension of the 
system Hilbert space and can be implemented with avail- 
able computer software.) Our second main result demon- 
strates that in certain cases finding active quantum error 
correcting codes reduces to the problem of finding pas- 
sive codes. Specifically, in the case of a noise operation 
that is unital but otherwise arbitrary, we identify a sub- 
set of the actively correctable subsystems, those that can 
be corrected by a single unitary operation, as precisely 
the noiseless subsystems for the composition of the op- 
eration with its dual. Hence the identification of such 
actively correctable subsystems may be accomplished via 
an application of the passive theory 0, |3] . 

Let us first discuss nomenclature. Given a quantum 
system S represented on a (finite-dimensional) Hilbert 
space Tt = Ti., we say a quantum system B is a subsys- 
tem of S if there is a representation of B on a Hilbert 
space Tt B such that 

H s = {H A ®H B )®K,, (1) 

where A is also a subsystem of S and K, = (H A ® Tt B ) ± . 
We are interested in cases dim7i s > 1, so at least one 
qubit can be encoded. We adopt the convention that p, 
a and r denote density operators and a superscript such 
as a B refers to the subsystem on which the operator is 
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defined. The set of operators on TL is denoted by C(TL). 

Linear maps on C(TL) can be regarded as operators 
acting on the space C{TL) with the Hilbert-Schmidt in- 
ner product (<t,t) = Tx(o"Tr). These maps are called 
superoperators to distinguish them from operators act- 
ing on TL. We use the terms channel and operation to 
mean a trace-preserving completely positive linear map 
£ : C(TL) — > C(TL') between Hilbert spaces H and TL'. 
For convenience we shall assume the noise operation £ 
maps TL to 7i, but our results are easily adapted to the 
general setting simply by replacing TL with TL' at the 
output stage. Such maps describe time evolution of open 
quantum systems (in the Schrodinger picture), and can 
always be represented in the operator-sum representation 
£( CT ) — J2aEaO~El by a set of "Kraus" operators {E a } 
on TL. The composition of two superoperators will be 
denoted by £ o T{a) = £(T~(a)). 

The dual or adjoint map £^ is the superoperator 
on C{TL) defined in the usual way by the equation 
Tr(£t(<7)r) = Tr(cr £(t)). A unitary operation U sat- 
isfies o hi = U o = id, where id is the identity 
superoperator. This is equivalent to the existence of a 
unitary U € £(7i) that implements the operation via the 
equation U(a) = UaW . Given a decomposition of TL 
as in Eq. (JJJ, Tab denotes the superoperator defined by 
"Pab(&) = Pab&Pab where Pab is the projector onto 
the subspace TL AB = TL A ® TL B , and let ids be the iden- 
tity superoperator on C(TL B ). If we are given maps £a, 
£b on the subsystems, as a notational convenience we 
write £a <8> £b both as a map on C(TL A <S> 7i B ) and for the 
natural extension of the map to C(TL). 

The definition from |3|, |4| of a correctable subsystem 
for a quantum operation £ is as follows. 

Definition: B is a correctable subsystem for an op- 
eration £ if there is a quantum operation TZ, which we 
call the correction operation, and a quantum operation 
T A ■ C(TL A ) -> such that 

Va A Vo- B : (fto£)(<7 A ®<7 B ) = JF a (ct a ) ® cr 8 , (2) 

or equivalently, such that 

TZo£oVab = Fa® ids- (3) 

It is worthwhile to define a variant of this notion. First 
observe that in general there will be many representations 
of B as a subsystem of S. That is, typically there will 
exist many subsystems C and Hilbert spaces KJ such that 

TL S = {TL C ®TL B )®K! . (4) 

Of course, C can be distinct from A, and can even be 
associated with a Hilbert space of different dimension. 
It is the consideration of different representations of a 
subsystem that leads to the following notion. 

Definition: B is a recoverable subsystem for an oper- 
ation £ if there is a quantum operation TZ, which we 



call the recovery operation, and a quantum operation 
T c \a '■ £{T~L A ) — ► C(TL°) where C is a subsystem of S 
that can be different from A such that 

V<t a V<j b : {no£)((T A ®o- B ) =F c \a(o A )®° B , (5) 

or equivalently, such that 

TZ o £ o Tab = Fc\a ® ids- (6) 

Note that a correctable subsystem is a special case of a 
recoverable subsystem where the subsystem C is equiva- 
lent to A. Furthermore, if B is a recoverable subsystem, 
then one can always compose the recovery operation TZ 
with an additional operation TZ' that maps all states of 
subsystem C to states of A such that the composition 
TZ' o TZ satisfies the definition of a correction operation. 
Consequently, B is a recoverable subsystem if and only 
if it is a correctable subsystem. (It is worth noting that 
the operation TZ' will typically correspond to "cooling" 
the system if the dimension of C is greater than that of 
A, however this can always be achieved by appending to 
the system an ancilla of dimension at most the square of 
that of the system and implementing a unitary operation 
on the composite HJ.) 

There is a significant difference, however, between im- 
plementing a correction operation and implementing a 
recovery operation. In both cases, quantum information 
stored in subsystem B can be brought back to B after 
a single iteration of the noise operation £, but only if 
one implements a correction operation will B be recov- 
erable from a subsequent iteration of the noise operation 
£ . The reason is that if one only implements a recovery 
operation, then the state is left in CB, while quantum 
information in B is only guaranteed to be recoverable if 
the initial state is in AB. 

Furthermore, when we constrain the nature of the re- 
covery and correction operations, the two notions no 
longer coincide. Specifically, we can introduce the fol- 
lowing two notions. 

Definition: B is a unitarily correctable subsystem for 
an operation £ if it satisfies the definition of a correctable 
subsystem and the correction operation TZ can be chosen 
to be a unitary U, so that 

Uo£oVab =T a ®\A b . (7) 

Definition: B is a unitarily recoverable subsystem for 
an operation £ if it satisfies the definition of a recoverable 
subsystem and the recovery operation TZ can be chosen 
to be a unitary hi, so that 

U o £ o Vab = Fc\A ® ids- (8) 

Note that by operating with W from the left on 
Eqs. J7J) and (JSJ), we obtain an expression for the ac- 
tion of the noise operation on AB. Heuristically, being 
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unitarily correctable or recoverable ensures that although 
quantum information stored in the subsystem B may be 
moved to another subsystem, it remains coherent. In the 
general case of possibly distinct input and output Hilbcrt 
spaces, we note that the relevant notion to consider is 
that of isometric recovery. 

Although a unitarily correctable subsystem is unitarily 
recoverable, the converse is not true unless C satisfies 
dimC = dim A because only in this case can one find a 
unitary operation TV that maps all states of subsystem 
C to states of A. Thus all we can say in this case is 
that a unitarily correctable subsystem is equivalent to 
a unitarily recoverable subsystem for which dim C is at 
most dim A. 

One might expect that not every correctable subsys- 
tem is unitarily correctable, and this is indeed the case. 
However, every correctable subsystem is unitarily recov- 
erable. The converse of this implication is also true, and 
together these facts constitute our first main result. 

Theorem 1 Let £ be a quantum operation. The follow- 
ing conditions are equivalent: 

(i) B is a correctable subsystem for £ . 

(ii) B is a unitarily recoverable subsystem for £ . 

The proof is presented in Appendix A. Note that it is 
a constructive proof in the sense that if an error model 
£ and a correctable subsystem B are known, then the 
unitary operation which yields Eq. JHJ may be explicitly 
obtained. Moreover, as discussed above, a correction op- 
eration can be obtained simply by composing U with an 
additional operation TV that maps states of C to states 
of A such that 71 = 11' oU satisfies Eq. ©. 

There are two special cases of this theorem that are 
worthy of note. First, consider the case where B is a 
correctable subsystem for £ and Ti is one-dimensional, 
so that Tt A <g> 7~t B is a correctable subspace for £. In this 
case, Tc\a is simply a density operator wc, and Eq. (jHJ) 
becomes Vct b : U o £{a B ) = lo c ® a B . Theorem Q] 
therefore establishes that a correctable subspace for £ is 
mapped by £ to a subsystem, modulo an overall unitary. 
This special case was proven by Nayak and Sen ■ That 
the subspace gets mapped to a subsystem may be viewed 
from our general perspective as simply a representation 
change of the subsystem B. 

The second special case of the theorem established pre- 
viously is the case H = H A ®H B , or equivalents K. = {0} 
in Eq. In other words, B is associated with a factor 
space of TL rather than a factor space of a subspace of 
H. In this case, any subsystem C in Eq. © must have 
the same dimension as the subsystem A. Consequently, 
there is no distinction between unitarily recoverable and 
unitarily correctable in this case. Dropping the projec- 
tion operation Vab, which is the identity operator in this 
case, Eq. © becomes U o £ = Ta ® ids or equivalently, 
£ = W o {Ta ® ids) • Theorem ^ therefore establishes 



that if a factor space B is a correctable subsystem for £, 
then £ simply maps B unitarily to another factor space. 
This case was proven by Nielsen and Poulin • 

Theorem ^ allows us to immediately make a nice con- 
nection with the testable conditions for operator quan- 
tum error correction. Taking the adjoint of Eq. (JSJ, we 
have 

Tab °£ f oW* =Tl lA ®id B . (9) 

Composing Eqs. (|7I9|I . noting that W oU — id and defin- 
ing the map Ga- B{Ua) -* B{Ha) by Qa = ^ C \ A °^c\a, 
we obtain 

Testable condition for B to be a correctable 
subsystem for £ : There exists a positive su- 
peroperator Ga on A such that 

V AB o£^ o£oVab = Ga®'^b. (10) 

This is simply the superoperator form of the standard 
testable condition for operator quantum error correction 
El E3- This is seen by noting that two completely 
positive maps are equal if and only if the Kraus oper- 
ators of one are related by a unitary remixing of the 
Kraus operators of the other |2JJ . If the Kraus operators 
of £ are denoted by {E a } and those of Qa are denoted 
{G c }, then the Kraus operators of Tab ° £' ° £ ° Vab are 
{PABElE a P AB } a ,b and those of GA®id B are {G e ®I B } c , 
so that Eq. I|1U|) is equivalent to 

P AB E\E a P AB = F ab ® I B (11) 

where F a = ^2 c u aiC G c is a unitary remixing of the G c . 
This is the testable condition discovered in 0, • It fol- 
lows from the results of [23| that not only can one derive 
Eq. I |10|l from the fact that B is a correctable subsys- 
tem for £, as we have done, but the opposite implication 
holds as well. (It would be interesting to have a version 
of this direction of the implication that is native to the 
superoperator formalism.) 

If A is one-dimensional, B is a subspace of S and the 
positive superoperator Ga is simply a positive scalar 7 2 , 
so that we have V B o £^ o £ oV B = "f 2 V B . This is the su- 
peroperator form of the testable condition for correctable 
subspaces. The condition was presented in this form in 
0. 

Finding correctable subsystems (or subspaces) for ar- 
bitrary noise models in full generality appears to be an 
intractable problem. Nevertheless, it has been possible 
to do so in several special cases, and this has been one 
of the successes of quantum error correction. The case 
of passive error correction is one such instance. As an 
application of Theorem ^ together with further analysis, 
we next show how recently developed techniques for pas- 
sive error correction 0, |2fJ can be extended to a special 
case of active error correction. The notion of a noiseless 
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(or decoherence-free) subsystem is fundamental to pas- 
sive error correction and so we begin with its definition. 

Definition: B is a noiseless subsystem for an oper- 
ation £ if there is an operation Qa ■ C(H A ) — > C(H A ) 
such that 

Va A Va B : £ (a A ® a B ) = Qa{<?a) ® ffB. (12) 
or equivalently, such that 

5 oP AS =0 A ®id B . (13) 

We will now focus on united quantum operations. Uni- 
tal operations are the quantum analogue of classical "bis- 
tochastic" maps (defined in Appendix B) and are ubiq- 
uitous in quantum information and computation A 
unital operation £ is one which takes the identity opera- 
tor to itself, £{I) = I. The dual of a unital map is neces- 
sarily trace-preserving because Ti(£^(a)) = Tr(£(I)cr) = 
Tr(<r). Thus, if £ is a unital channel, then f ' is a can- 
didate for a correction operation. Moreover, the dual of 
a trace-preserving map is necessarily unital because if £ 
is trace-preserving then Tr(£T (I)cr) = Tr(£ (a)) = Tr(cr) 
for all a and consequently £ = I. Thus, if the super- 
operator £ is both unital and trace-preserving, then so is 
£' and consequently so is £ * o 5. 

It is easy to see that if B is a noiseless subsystem for 
£ ' o 5, so that f'ofo "Pas = £m <8> ids, then B is a cor- 
rectable subsystem for 5 because £^ is trace-preserving 
and thus constitutes a correction operation that satisfies 
the definition of Eq. . (The correction in this case sim- 
ply involves "reversing" the noise on the system.) How- 
ever, one can say more than this, namely, that B is uni- 
tarily correctable for £ . Together with its converse, this 
fact constitutes our second main result. 

Theorem 2 Let £ be a unital quantum operation. The 
following are equivalent: 

(i) B is a unitarily correctable subsystem for £ . 

(ii) B is a noiseless subsystem for £ t o £ . 

The proof relies on a number of ancillary results for 
unital operations presented in Appendix B. The following 
lemma, which applies to arbitrary positivity-preserving 
superoperators, does not rely on these results and hence 
we prove it here. 

Lemma 3 Let O be a positivity-preserving superopera- 
tor. Let P be a projector and let V(-) = P(-)P be the 
associated projective superoperator. Then 

VoOoV = OoV (14) 

if and only if 

supp(0(P)) C supp(P) , (15) 
where supp(er) denotes the support of the operator a. 



Proof. The superoperator identity Eq. (|14|l implies 
PO(P)P = 0(P) which is equivalent to Eq. |JT5J}. Con- 
versely, suppose Eq. Ijl5|l holds and let p be an arbitrary 
density operator. Noting that P — V{p) is a positive 
operator, it follows from the fact that O is linear and 
positivity-preserving that O(P) -Co is also pos- 
itive and consequently supp(C °V(p)) C supp(0(P)). 
Given Eq. JT2J), we have supp(C o V{p)) C supp(P). Fi- 
nally, given that p is arbitrary, Eq. (|14|) follows. QED 

The following lemma, the proof of which can be found 
in Appendix B, is also central to the proof of the theorem. 

Lemma 4 Let £ be a unital operation, and let B be a 
correctable subsystem for £ . The following are equivalent: 

(i) supp(£ t o £{P AB )) C supp(Pab) 

(ii) rank(£ (Pab)) = rank(PAs) 

We now provide the proof of Theorem [3 
Proof of Theorem We first establish the implication 
(i) =>• (ii). Given that B is a correctable subsystem, 
it follows that the testable condition of Eq. (|10|l holds, 
that is, Tab o £^ o £ o Tab = Qa® ids (and Theorem^] 
provides a simple way of seeing this). At this stage, it is 
clear that if not for the leading Tab in this expression, 
B would satisfy the definition of a noiseless subsystem 
for £^ o £ provided in Eq. (|T3|> . namely, (£t o £) o Tab = 
Qa <8> ids • By Lemma |3| one sees that the leading Tab 
can be dropped if supp(£t o £(Pab)) Q supp(Pab), so it 
remains only to show that the latter condition is satisfied 
if B is unitarily correctable. 

Since B is a unitarily correctable subsystem for 
£, there exists a unitary operation IA such that 
U o £ o Tab — 3~a ® ids, from which one 
finds that £{P AB ) = UHT A (I A ) ® Ib). It fol- 
lows that rank(£(P4s)) = rank(Jp4(l4))rank(/B) < 
rank(/ J 4)rank(/s) = rank(P4s), where we have made 
use of the fact that Pab — Ia® Ib- But £ is unital, 
so by Corollary (provided in Appendix B), we have 
ra,nk(£ (Pab)) > rank(PAs). One concludes that there 
must be equality between the two quantities, that is, 
rank(£(P4s)) = rank(P4s). By Lemma 0J it follows 
that suppO? 1 o £(P A b)) Q supp(P4 B ). 

We now prove (ii) =4> (i). From the fact that 
B is a noiseless subsystem, we infer that supp^ o 
£(Pab)) Q supp(Pab). By Lemma it follows that 
rank(£(PAs)) = rank(PAs). Now, given that B is a 
correctable subsystem for £ (as established in the dis- 
cussion above Theorem EJ, by Theorem ^ it is also a 
unitarily recoverable subsystem for £, so that there ex- 
ists a unitary operation U such that U o £ o Tab — 
T C \a ® ids- Given that Pab — I a® Ib, it follows that 
U o £(P A b) = 3~ C \a(Ia) ® Ib, and consequently that 
rank(£(P4s)) = rank(J r C | j 4(/ y i))rank(/B). But given 
that rank(£(P4s)) = rank(P J 4s) = rank(l4)rank(Is), 
we conclude that rank^c^T/i)) = rank(Pi). This im- 
plies that there exists a unitary V such that V o (!Fc\a ® 
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ids) = JF A (g)id_B for some operation Ta on subsystem A. 
Consequently, VoUoEoVab = VoTc\a®'^b = Ta®'^b, 
which by Eq. Q implies that the unitary operation VoU 
is a correction operation for £, and thus that B is uni- 
tarily correctable. QED 

By combining Theorem [S] with the recentl y d eveloped 
structure theory for noiseless subsystems [il l2flj. we ob- 
tain a method for finding the unitarily correctable subsys- 
tems for any unital operation £. Specifically, the noiseless 
subsystems of £^ o £ are obtained from the fixed point 
set Fix(£T o £) — {a : (£* o £)(a) — a}, in the following 
way. This set is a f-algebra [2j| and the representation 
theory for such algebras induces a Hilbert space decom- 
position TL = ®k(H Ak <8> Tt Bk ) in which B is a noiseless 
subsystem for £^ o £ if and only if TL B C 7i Bfc for some 

[2(| • (See 0, |2(| for further discussions and analy- 
sis.) The fixed point set of an operation £ is simply the 
eigenvalue-1 operator eigenspace of the superoperator £ , 
which is straightforward to determine. 

For an arbitrary unital noise operation £, not every 
correctable subsystem is unitarily correctable, and con- 
sequently the noiseless subsystems of £> o £ do not in 
general capture all correctable codes for a typical unital 
channel £ . 

The generic two-qubit "binary unitary channels" pro- 
vide a class of operations that illustrate this point [27j. 
As an example, let U be a unitary on TL = C 4 = C 2 (g> C 2 
with distinct eigenvalues Xj = cxp(i#j) ordered so that 
< 6»i < 6 2 < 9 3 < 04 < 2vr. Let j = 1,2,3,4, be 
corresponding eigenstates. Fix a probability < p < 1, 
and define a unital channel £ by £ (a) = pa+(l—p)UaU'. 
Then £ has two Kraus operators given by, up to normal- 
ization, £ = {/,[/}. The noiseless subsystems (actually 
subspaces in this case) for £T o£ come from the so-called 
noise commutant {Pi, PU, Wl, WU}' = {U, W}' = 
{[/}', which coincides with the fixed point set. But since 
U has distinct eigenvalues, this commutant is isomorphic 
as a f-algebra to the algebra C C C C, and hence 
can only be used to encode classical information. 

On the other hand, this channel has correctable qubit 
codes. For example, let A be the point of intersection 
of the line segments [Ai,Aa] and [A2,A4]. Let s and t 
be fixed probabilities such that A = sAi + (1 — s)A3 = 
t\ 2 + (1 — t)A 4 , and define states {|^>), \(f>)} by 

|V) = \th) + Vi^s |v 3 ) 

10) = Vt\i/J2) + VT~t\i; 4 ). 

The two-dimensional subspace spanned by and \tp) is 
a correctable code for £ . In particular, one can compute 
that PUP = \P where P = \tp)(ip\ + \(/))((f>\, and thus 
the error correction condition from 0] is satisfied for £ 
on the subspace C = span{|?/>), \ <p)}. 

There are many unital noise operations, however, for 
which the composition of this map with its dual does have 



noiseless subsystems. As a simple example, consider the 
swap operation \ip) ® \<t>) >— » \4>) ® 1^) on a composite 
quantum system H — H A (8 TZ A made up of a subsystem 
TL A and a replication 1Z A = TL A . It is clear that both 
the subsystem 7i A and its copy can be returned to their 
initial locations by simply applying the swap operation 
again (which is equal to its dual). 

Of course, one could note that the swap operation itself 
has a noiseless subsystem of the same size; namely the 
symmetric space \ip)®\ip). But it is easy to find examples 
of operations with no noiseless subsystem, for which the 
composition map has a non-trivial noiseless subsystem. 

To this end, consider a two-qubit system exposed to de- 
coupled phase flips. The associated error model satisfies 
£{p) = pZ\pZ\ + (1 — p)Z 2 pZ 2 for some fixed probability 
< p < 1 and Zi = Z <g> il 2 , Z 2 = 11 2 ® Z. In this case 
there is no noiseless subsystem (or subspace) for £ . This 
follows from the fact that the noise commutant {Z 1: Z 2 }' 
is isomorphic to the algebra C C C C. Thus, only 
classical information can be safely sent through the chan- 
nel unscathed. 

However, the operators supported on the subspace 
spanned by \0l) = |00) and |1l) = |11) form a noise- 
less subspace for £^ o £. Indeed, the set of opera- 
tors a = a|00)(00| + 6|00)(11| + c|ll)(00| + d|ll)(ll| 
form a subalgebra of the commutant Fix(£T o £) = 
{Z{Z 2 ,z{Zi,zlZiZlz 2 }' = {ZiZ 2 }'. The unitary cor- 
rection operation guaranteed by Theorem [21 in this case 
happens to be the controlled phase flip operation U — 
|00)(00| + |01)(01| + |10)<10|-|11)<11|. 

Conclusion. — We showed that every correctable sub- 
system for an arbitrary quantum operation is a unitarily 
recoverable subsystem for the operation. Thus, the ef- 
fect of the operation on the subsystem can be reversed 
by a single unitary operation, up to a change in the rep- 
resentation of the subsystem. Our proof was construc- 
tive in nature, showing explicitly how the representation 
and unitary may be obtained, and hence a correction op- 
eration, if the operation and correctable subsystem are 
known. We also suggested that this result contributes to 
the rigorous formulation of the subsystem principle for 
quantum operations. 

We showed that the unitarily correctable subsystems 
for unital quantum operations are precisely the noiseless 
subsystems for the operation followed by its dual. We 
indicated how such subsystems can be practically com- 
puted and discussed some simple examples. The possibil- 
ity of extending this result to the case of nonunital quan- 
tum operations remains a problem for future research. 
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APPENDIX A: PROOF OF THEOREM [T] 

Proof of Theorem^ The implication (ii) => (i) follows 
immediately from the fact that a unitarily recoverable 
subsystem is a special case of a recoverable subsystem, 
and as discussed below Eq. JSJ being recoverable implies 
being correctable. 

For (i) (ii), let £ = {E a } be a Kraus operator 
representation for £ and assume B is correctable for £. 

We proceed as follows. First we shall construct an op- 
eration Q = {G a } such that the operators {G a P AB } have 
mutually orthogonal ranges, that is P AB G\G b P AB = 
for all a ^ b, and Q{I A ® a B ) = £{I A <g) a B ) for all 
a B . Then we shall find a subsystem C such that CB 
can be identified with a subspace of TL, an operation 
Tc\a '■ C(TL A ) — > C(H C ), and a unitary operation V 
on C(7i), such that 

g{a A ®a B ) = V{T c]A (a A )®a B ) V<j a Va B . (16) 

Then we will have 

V^o£(I A ®o- B ) = T c]A (I A )®a B V<t b , (17) 

and we can use a customary positivity-cum-linearity ar- 
gument (c.f. Lemma 2.3 4]) to show that Eq. © holds 
for all a A . 

We begin by noting that the testable conditions from 
0, 0, H3| give us operators F ab on 7i A such that 

P AB EiE b P AB =F ab ®I B Va,b. (18) 

Observe that the operator block matrix F = (Fat,) is 
positive since (I m ®P AB )E' s E(I m ®P AB ) = F®I B , where 
the row matrix E = [E ai E a2 ■ ■ ■], the number of E a is to, 
and I m is the identity operator on TO-dimensional Hilbcrt 
space. 

Next let U be a unitary such that U FW = D is diag- 
onal, and let U = (U a b) and D — (D a b) be the associated 
block decompositions. Then 

J2 U ^F c dU t bd = d ab D aa Va,6, (19) 

C,rf 

UlU cb = S ab I A Vo,6. (20) 

c 

Define a superoperator Q — {G a } where for all a, 
G a = J2 E ^ U lb ® i b)Pab + E a Pi B . 

b 

Let X ab = E b {Ul b <g> I B )P AB . Then by Eqs. CEEHj , one 
can verify for all a, b 

PABG\c b p AB = Xl c X bd 

c,d 

= (j2u ac F cd ul) ® I B 

c,d 

= D ab ®I B , 
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and D ao — for all a ^ b. Moreover, Eq. 1|20[1 yields for 
all a B 

G{I A ®a B ) = Y,Ga(lA®v B )Gl 

a 

= Y. X *b{lA®o B )Xl 

a,b,c 

= Y. E ^iX. U lbUac)®OB)El 

b,c a 

b 

= £{I A ®a B ). 

To simplify notation, let dim B = n and fix an or- 
thonormal basis {\ipk)}k=i f° r B. Let m be the car- 
dinality of the set {a : D aa ^ 0}, and for all a let 
r a = rank y/D aa = rankP aa . 

By the polar decomposition applied to each G a P A s, 
and the fact that these operators have mutually orthogo- 
nal ranges, there are partial isometries V a with mutually 
orthogonal ranges for distinct a such that 

G a P AB = V a \j PABGaGaPAB = V a {^D~ ® I B ). (21) 

It follows that we can find a decomposition of each V a 
of the form 

Tl r a 

^ = EEK (a) )K a) )^i^ (a) i' ( 22 ) 

k=l 1 = 1 

where {|</>| a ' ) )}[=i is an orthonormal basis for 
range \J D aa = range D aa for each a. The vectors 
I^W^gW^ _ V a \(j>i)\ipk) are given by the tensor 
product on the final space V a V a ^TL induced by the tensor 
product on the initial space VJV a TL and the isometric 
action of V a : VjV a TL -> V a V^H. 

Now let {|w a )}aLi be an orthonormal set of vec- 
tors in an m-dimensional Hilbert space. The vectors 
{\(f>l}\w a )\ipk}} a ,k,l are orthonormal since the \w a ) are, 

and because, for a fixed a, the vectors {|</ , |°'')|V-'fc)}fc,i form 
an orthonormal set. Moreover, the orthogonality of the 
ranges of the operators G a P A s ensures the dimension 
of TL is at least X)!ILi nr a- Thus, we may identify the 
Hilbert space spanned by the vectors {|0[ a ') \w a )\i^).)}a k,l 
with a subspace of TL. 

Next define V to be any unitary extension of the fol- 
lowing partial isometry to all of TL: 

V : \4 a) )\v) a )\fl> k )^\ir^)\6^). (23) 

The image vectors {\tt\ \)\0j^)}a t k,l form an orthonormal 
set since the V a are partial isometries with mutually or- 
thogonal ranges, and because these vectors are orthonor- 
mal for a fixed a. Further define Tq\a = {\/D aa ® \w a }} 



as a channel on C(TL ), and so for all a A , 

^c|a(0 = E(v^* a a/^)®K)K|. (24) 

a 

Finally, it can be verified by direct computation that 

g{a A ®u B ) = V(T c \ A {a A )®a B ) Va A Va B , (25) 
and this completes the proof. QED 

APPENDIX B: THEOREM [2] 

In this appendix we establish results used in the proof 
of Theorem [2] 

We begin by recalling a result for classical maps. A 
vector of probabilities p of dimension n is said to be 
majorised by a vector q, denoted p -< q, if for each k in 
the range 1 to n, Ylj=iPj — Ylj=i Iji wrtn equality for 
k = n, where the J, indicates that the probabilities are to 
be taken in nonincreasing order. A bistochastic map is a 
matrix A satisfying J2 k Akk' — 1 and J\' ^fcfc' = 1. The 
Hardy-Littlewood-Polya theorem [2^, l30| states that the 
output of a bistochastic map is majorised by the input, 
that is, if p = Aq for bistochastic A, then p -< q. 

A similar result holds for unital operations. It is a 
theorem due to Uhlmann [3ll l32| . We include a short 
proof for completeness. 

Lemma 5 If p = £(o~) for a unital channel £ , then the 
ordered spectrum r of p is majorised by the ordered spec- 
trum s of a; that is, 

p = £(o~) with £ unital implies r -< s. (26) 

Proof. Let pi~ and \e k ) (respectively, and \fk)) de- 
note the fcth eigenvalue and normalized eigenvector of 
p (respectively, a). Clearly, q k = J2k> D kk'Pk> where 
D k k' = Tr(|/ fe ) {f k \£(\e k ') {e k '\)) ■ From the fact that £ is 
trace-preserving one infers that Dkk' = 1 while from 
the fact that £ is unital one infers that J2 k ' Dkk' = L 
It follows that D = (Dkk 1 ) is bistochastic. Thus, by the 
Hardy-Littlewood-Polya theorem, we have r -< s. QED 

Heuristically, this says that unital channels can only 
increase the impurity or "mixedness" of quantum states. 

Corollary 6 For a unital operation £ acting on a pro- 
jector P, we have rank(£(P)) > rank(P). 

This is a simple consequence of the majorisation rela- 
tion given by Eq. (12611 . 

Corollary 7 Let £ be a unital operation and P a projec- 
tion. 7/rank(£(P)) = rank(P). that is, if the inequality 
of Corollary\Qis saturated, then £(P) is a projection. 
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Proof. A projection P has a uniform spectrum, and the 
only spectrum that is majorised by the uniform spectrum 
while having the same number of non-zero elements is 
the uniform spectrum. Given that £ is trace-preserving, 
it follows that £(P) is a projection. QED 

Lemma 8 For a unital operation £ acting on a projector 
P, we have £{P) = P if and only iftf(P) = P. 

Proof. By the fixed point theorem for unital opera- 
tions |2||, P is in the fixed point set of £ if and only if 
it commutes with all Kraus operators of £. The lemma 
then follows from the fact that the Kraus operators of 
are the adjoints of those of £ and that both maps are 
unital. QED 

We are now in a position to prove Lemma 01 
Proof of Lemma^ For (ii) (i), we assume that 
rank(£ (Pab)) — rank(PAs). It follows from the fact that 
£ is unital and Corollary[7] that £{Pab) is a projection. 
A projection of the same rank as Pab is obtained by a 
unitary operation IA acting on Pab, so that £(Pab) = 
U(Pab)- We infer that U ] o£(P AB ) = P AB . By LemmaEl 



it follows that £' oU(Pab) — Pab- Consequently we have 
£t o£(P AB ) = £' oU oW o£ (Pab) = Pab, and a fortiori 
condition (i) holds. 

For (i) (ii), we assume that supp(£^ o £(P A b)) Q 
supp(Pab)- It follows that rank(£t £(P A b)) < 
rank(PAs)- Thus, given that £ ' o £ is a unital opera- 
tion, by Corollary El we deduce that in fact rank(£t o 
£(Pab)) — rank(PAs)- By Corollary |7| it follows 
that £t o £(P AB ) is a projection. The only projec- 
tion with the same rank and support as Pab is Pab, 
therefore £ j; o £(P AB ) = Pab- Given that £t is uni- 
tal, it follows from CorollarylfJlthat rank(£t o£(P AB )) > 
ra,nk(£(P AB )). Given that £ t o £(P AB ) — Pab, we de- 
duce that rank(PAs) > rank(£ (Pab))- However, from 
Corollary El and the unitalness of £, we infer the opposite 
inequality rank(£(P J 4s)) > rank(P4s), so that we must 
have equality, ia,nk(£ (Pab)) — rank(PAs). QED 

Notice from this proof that the two conditions of 
Lemma 0] are equivalent to the seemingly stronger state- 
ment £t o £(P A b) = Pab- 



